ANNENBERG 


BROWN UNIVERSITY 


sXe WYAVCol a dlale] at-) ol-1am \ opal r-ts10) 


Gender Peer Effects in Post-Secondary 
Vocational Education 


Fernanda Ramirez-Espinoza 
Harvard University 


This paper presents evidence that women and men benefit from having a higher percentage of female peers in 
post-secondary vocational STEM programs. I use idiosyncratic variation in gender composition across 
cohorts within majors within branches (campuses) for identification. Having a higher percentage of female 
peers positively affects students in STEM majors, decreasing women's dropout rates and increasing GPA. The 
peer effect seems to be mediated by the gender of the instructors: as female students have fewer female 
instructors, the effect of having more female peers intensifies. For men, a higher percentage of female peers 
reduces dropouts and increases GPA to a lesser extent, suggesting that policies that increase the representation 
of women need not entail a trade-off for male STEM students. 


VERSION: October 2021 


Suggested citation: Ramirez-Espinoza, Fernanda . (2021). Gender Peer Effects in Post-Secondary Vocational Education. (EdWork- 
ingPaper: 21-480). Retrieved from Annenberg Institute at Brown University: https://doi.org/10.26300/s6r2-sw83 


Gender Peer Effects in Post-Secondary Vocational 


Education ! 
Fernanda Ramirez-Espinoza * 


September 29, 2021 


Abstract 


This paper presents evidence that women and men benefit from having a higher 
percentage of female peers in post-secondary vocational STEM programs. I use idiosyn- 
cratic variation in gender composition across cohorts within majors within branches 
(campuses) for identification. Having a higher percentage of female peers positively 
affects students in STEM majors, decreasing women’s dropout rates and increasing 
GPA. The peer effect seems to be mediated by the gender of the instructors: as female 
students have fewer female instructors, the effect of having more female peers intensi- 
fies. For men, a higher percentage of female peers reduces dropouts and increases GPA 
to a lesser extent, suggesting that policies that increase the representation of women 
need not entail a trade-off for male STEM students. 


Although educational attainment gaps have not only narrowed but have reversed in most 
high-income countries and Latin America (Goldin, 2002; Goldin, Katz, & Kuziemko, 2006; 
Duryea, Galiani, Nopo, & Piras, 2007), a high degree of occupational segregation remains: 
men and women are still concentrated in different occupations (Schneeweis & Zweimiiller, 
2012). This is an essential point in that gender wage differences are partly attributable to 
the subjects that men and women choose to study. Consequently, studying the mechanisms 
through which women persist (or desist) in high-paying paths like STEM fields is relevant 
as it helps us understand how to close this persistent gender inequity. 

In this paper, I hypothesize that female students’ educational outcomes will be positively 
related to having more female peers in their cohorts in STEM majors. My analysis is moti- 


vated by previous literature that suggests that peer effects exist and are particularly salient 
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in the STEM fields, where females have lower participation than in other fields (Hoxby, 
2000; Bostwick & Weinberg, 2018; Sacerdote, 2011). The identification strategy to test this 
hypothesis uses observational data from 96.627 students of the largest post-secondary vo- 
cational institution in Chile. The outcome studied is first-year dropout and Grade Point 
Average (GPA). 

To identify the gender peer effect, I estimate a linear model of educational outcomes on 
gender composition and gender of the student, including an interaction term that captures 
how female students differentially respond to gender composition. To control for unobserved 
characteristics of the students that might be related to dropout rates and gender composition, 
I rely on deviation from long-term trends in gender peer composition within major-by-branch 
across the years. For example, to calculate the effect, I compare different cohorts of Automo- 
tive Mechanic in the Valparaiso campus, where students were exposed to a slightly different 
percentage of females year to year. To do so, I include a major-by-branch fixed effect and 
time trend. Additionally, I include control variables such as age, diagnostic math test scores, 
education of the mother, working status, financial aid status, and shift (day/night). 

Many researchers and teachers have argued that peer composition is as important a 
determinant of student outcomes (Sacerdote, 2011). This issue seems to be particularly 
crucial for females: there is evidence that females respond more than males to peer influences, 
consistent with social psychology theories that peers affect female students more (Han & Li, 
2009). 

There is evidence that the gender composition of peers affects outcomes and that these 
effects are different for boys and girls (Busso & Frisancho, 2021; Mouganie & Wang, 2020; 
Zolitz & Feld, 2020). Lavy and Schlosser (2011) find large positive effects from the percent 
of girls within a classroom, and they also interpret these effects as working through more 
than merely increasing peer average test scores. Along the same lines, (Hoxby, 2000) finds 
modestly large effects of peer background on own test scores, using idiosyncratic gender 


variation. Paredes (2018) finds that single-sex classrooms reduce the math gender gap by 


more than half in the Chilean context. She finds that this effect is driven by the gender 
composition of the classroom itself. 

Although the gender peer effect literature in primary and secondary education is robust, 
post-secondary education remains a less-explored area. Some papers that study how culture 
may be connected to female underrepresentation have been published recently, like those 
authored by Lundberg (2017) and Wu (2017). Likewise, there are some studies on STEM 
graduate program admissions and persistence. For instance, Bostwick and Weinberg (2018) 
use a difference-in-difference approach and find that an increase in the percentage of female 
students differentially increases the probability of on-time graduation for women. This pa- 
per contributes to the emergent body of evidence of gender peer effects in post-secondary 
education. 

Furthermore, studies on gender peer effects in vocational education are almost non- 
existent, although it comprises a significant part of educational systems worldwide. In some 
developed countries, one-quarter of cohorts pursue professional programs. In the United 
States, certificate graduation rates are burgeoning — tripling in recent years (Skills beyond 
school: synthesis report, 2014). This study is situated in Latin America, where the post- 
secondary education sector is also growing fast. This trend is observed in countries with the 
highest secondary education completion rates, such as Colombia, Mexico, Brazil, Chile, and 
Peru. As other Latin American countries raise their secondary education completion rates, 
this paper sheds light on how gender peer effects play out in this specific context. 

The results I present in this paper suggest that a 10% increase in the percentage of female 
peers within major-by-branch units, close to the mean idiosyncratic variation in the data, is 
associated with a reduction of 9.5% (1.6 percentage points) in female students’ dropout rate 
and a 0.05 standard deviations increase in GPA. This result supports the hypothesis that 
female students’ educational outcomes are positively related to having more female peers in 
their cohorts in STEM majors. For males, this relationship is of a smaller magnitude for 


both dropout and GPA and significant for the case of GPA, suggesting that men in STEM 


programs also benefit from having more female peers. These results are robust to different 
controls specifications and hold if estimated with a subsample of only majors with many or 
few students. 

Some may argue that the problem is not some particular gender dynamic in STEM 
majors but male-concentrated majors. If this were the case, we would observe a similar 
effect for male-concentrated STEM and Non-STEM majors. Nevertheless, when the effect 
is calculated for non-STEM majors that are male concentrated, we do not find a significant 
effect. 

The causal identification strategy relies on the assumption that the variation in gender 
peer composition within major-by-branches is as good as random. ‘To provide evidence 
supporting this idea, I estimate an autoregressive model of gender peer composition with 
major-by-branch and year-fixed effects. In this context, peer composition in the previous year 
is not significantly correlated to gender peer composition in the current year, supporting the 
idea that within major-by-branch, the change of percentage of female peers is idiosyncratic. 

The paper is organized as follows. Section 1 summarizes the key features of the appli- 
cation process in the vocational education institution that will be studied, describes data 
sources, reports summary statistics, and tests for balance on the treatment variable and the 
plausibility of the identifying assumption. In Section 2, I describe the identification strategy 
used in this study. I present the main results in Section 3 and provide a preliminary mech- 
anism exploration in Section 4. I conclude in Section 5 by interpreting my findings in the 


context of gender peer effects. 


1 Context and Data 


1.1 Institutional Context 


Chile’s education system comprises four levels of education: pre-school, primary education, 


secondary education, and tertiary education (also known as higher education). Three types 
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of institutions can provide higher education: Universities, Professional Institutes (IPs), and 
Centers for Technical Training (CFTs). The most significant difference between universities 
and the vocational sector is the type and length of training they provide. Universities focus 
on formal academic training, while [Ps and CFTs focus on developing practical work skills. 
This is reflected in the length of courses — the average minimum length of university degrees 
for incoming students is about nine semesters (Arango, Evans, & Quadri, 2016). Students 
choose a major when they are admitted, and most of their classes are with peers of the same 
or similar majors. 

This paper uses information from DUOC UC °, the largest vocational higher education 
institution in Chile. They tend to 19.3% of all vocational higher education students in the 
country, offering 75 majors in 9 areas of health, tourism, construction, and management. In 
2018, they had 102,817 enrolled students in their 15 branches, present in three regions of the 
country. Twelve of these branches are located in the Greater Santiago Area (Figure 1), which 
covers an area similar to New York City and has 6,257,516 inhabitants. Given that some 
majors are taught in more than one branch, there are 323 Major-by-Branch combinations. 

Although academic requirements to study in DUOC UC are not strict, the institution 
holds the highest level of accreditation in the system — an honor shared only with three of 
the best universities in the system. According to this metric, Duoc UC is the highest quality 


vocational higher education institution in the country. 


1.2 The enrollment process 


DUOC UC has a rolling admissions process that is first-come, first-served.” Besides some 
minimal academic requirements, there are no requisites for enrollment besides proof of pay- 
ment (first come-first served). Therefore, neither the students nor the institution can forecast 
the percentage of female peers accurately they will have at the moment of enrollment. In 


theory, only the last student who enrolls could know how many female peers she/he would 


3For more information, visit: http: //www.duoc.cl/international-affairs 
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Figure 1: Map of DUOC branches in the Greater Santiago area 


be exposed to, but the characteristics of who are already enrolled are not public, so the 
student does not have this information available when she/he makes the decision. Therefore, 
individuals students cannot know exactly how many female peers they will have. 

Another essential institutional detail is how higher education studies are taught in the 
Chilean context. For the higher education system in general and DUOC UC in particular, 
students enroll directly into specific majors. This has two significant implications: first, if 
they wish to change majors, they must re-enroll to the new major and start this second 
major from scratch. The second, and most important in the context of the study, is that 
they always study with the same group of peers: those from their same year-of-entry, same 
major and same branch. Therefore, the concentration of female peers is the same for all 
students in that major in a particular branch, and the opportunity for inter-major gender 
peer influence is less likely. 

As it was previously stated, DUOC UC has 15 branches in different parts of the country. 


Each of them offers a set of majors (e.g., Electrical Technician, Dental Technician, Gastron- 


omy), and most of them overlap. For this study, I consider the relevant peer group as those 
studying the same major at the same branch in the same year of entry. Therefore, a first- 
year Electrical Technician student in branch A has a different peer group than a first-year 


Electrical Technician student in branch B. 


1.3. Data 


I use data from all enrolled students in Duoc UC from 2014 to 2018. The dataset includes all 
individuals that studied one of the 76 majors continuously offered between 2014 and 2018 in 
its 15 branches. This dataset contains information on 104,146 students, with characteristics 
like gender, age, and mother’s education. It also includes information on their working 
status, the scores on a math diagnostic test all students take before the beginning of their 
first academic year, if they attend school during the day or at night (night shift), gender 
of the instructor, dropout, and GPA. All this information allows testing for the effects and 
mechanisms described in Section 2. 

The measurement of educational outcomes for this dataset is students’ dropout during 
her first year of studies at Duoc UC and students’ GPA. From the data, we get that 16.4% 
of students dropout in the first year. As shown in Table 1, STEM majors have a higher 
dropout rate than non-STEM majors (18.1% vs. 15.2%). STEM majors are markedly male 
concentrated: on average, they have 12.5% of females, whereas non-STEM majors have 60%. 
Compared to the mean, differences in the other covariates are small in magnitude, showing 
that the differences are not practically meaningful. 

The “treatment” in this setting is each student’s percentage of female peers in their first 
semester. This is calculated taking the number of female peers over the number of total 


peers per major-by-branch on their first semester: 


Lily; Female; 
PercentageFemale = ~ eae (1) 
Ntk — 


Table 1: Summary Statistics for STEM and Non-STEM students 


STEM Non-STEM _ Difference 
Dropout 0.181 0.152 Q.029°*% 
(0.002) 
GPA (in SD) 0.11 0.148 -0.258**# 
(0.006) 
Percentage Female 0.125 0.6 -0.475*** 
(0.001) 
Age 21.243 921.632 -0.389*** 
(0.113) 
Diagnostic score 48.819 45.936 2.883*** 
(0.121) 
Mothers’ education 0.632 0.639 -0.007** 
(0.003) 
Works 0.572 0.558 02444" 
(0.003) 
Has financial aid 0.661 0.666 -0.005 
(0.003) 
Night Shift 0.366 0.254 Caria 
(0.003) 
N 44,234 59,912 104,146 


Notes: This table presents means and mean difference 


Standard errors are reported in parenthesis (*p<0.1;**p<0.05; 


**<0.01) 
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Figure 2: Within Major-by-Branch variation in percentage of female peers (Valparaiso 
branch) 


Where n;, is student i’s number of peers n in the major-by-branch k in cohort t. 

The identifying variation used in this paper is the variation of the average percentage 
of female peers students have within the 323 major-by-branches. Figure 2 illustrates the 
variation of treatment within major for the one particular branch (located in Valparaiso). 
The year-to-year variation plotted in Figure 2 is the identifying variation in this context. 


The variation for all the institution branches can be seen in Figure A. 


2 Causal Identification Strategy 


In this context, the treatment is the percentage of female peers a student has in their program, 
in their branch, in their entrance cohort. Through a linear model that includes year-fixed 
effects, major-by-branch fixed effects, and time trends, I use year-to-year deviations from 
long-term trends in the percentage of female peers within Major-by-Branch to estimate the 
effect of interest. 


As the treatment was not randomized, there is the potential for bias in estimating the 


treatment effect. If the change of percentage of women were somehow correlated with unob- 
servables correlated with the outcome, the estimate of treatment effect would be biased. An 
example of an unobservable that could threaten the identification strategy is a program di- 
rector concerned about gender imbalance and makes efforts to attract more females to STEM 
programs and implement initiatives geared to eliminate the obstacles that make women drop 
out or get lower grades. 

In this case, the identifying assumption is that the treatment was assigned as if randomly 


conditional on the controls. If this is satisfied, the average treatment effect will be given by: 


ATE = E [UM (Yi (a fem+1%)/Xitk) = L (Yir(% fem) /Xitk) | (2) 


Where ¥it(¢% fem) is the percentage of female peers student i has in her cohort t (treatment), 
and Xj, are observable characteristics of student i of cohort t in major-by-branch k. Because 
the treatment Yi(q fem) 18 continuous, instead of using the traditional binary treatment of 
Yi(%@fem) being equal to 0 or 1, I use the continuous version of finding the difference in 
outcome between Yica fem+1%) ad Yic% fem) Conditioned in all other observables. 

To estimate the effect, I use a linear regression model. Because I am interested in knowing 
the effect of the percentage of female peers and how this effect is different for males and 
females, I will use an interaction term* between gender (male) and percentage of female 
peers. The linear regression model I will use to estimate the effect of the percentage of 


female peers on females and males is: 


“It is crucial to address the fact that when I include an interaction term with gender in the potential 
outcomes framework, there is an implication that both percentage of female peers and gender can be manip- 
ulated, analogous to treatment in a randomized experiment. If gender is not manipulable, then there is the 
danger of post-treatment bias stemming from the fact that almost all variables on which I will condition are 
determined after an individual’s conception (Greiner & Rubin, 2011). A shift in focus from actual traits to 
perceptions can address this problem, as suggested by Greiner and Rubin (2011). Here, the treatment is not 
a gender switch from female to male, but a change in the perceptions about students’ performance attached 
to their gender. As in a potential outcomes framework, we usually think of a state of the world we wish to 
achieve and an intervention that will get us closer to it. In this context, the goal would be that perceptions 
on performance are not related to a person’s gender, and the mechanism proposed to achieve such a goal is 
changes in the percentage of female peers the student is exposed to. 
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Yirx = Bo + Bi % female, + BoMale jz. + B3% female; * Male ig + BaX ink + Ye + Oe + Vavearst + Eitk 
(3) 

Yjx is the outcome of interest (dropout or standardized GPA) for student i in cohort t 

in major-by-branch k, % female, is the percentage of female peers in major-by-branch k in 
cohort t, Malej, is a dummy that takes value of 1 if the student i in cohort t in major-by- 
branch k is male, Xj, is a vector of student’s covariates, Y; are years fixed effects, 6, are 
major-by-branch fixed effects and wy; is a set of major-by-branch-specific linear time trends. 
The coefficient of interest is Bj, that indicates how percentage of female peers is related to 
dropout for females. If the assumptions hold, then B; would identify the causal effect of 


percentage of female peers on educational outcomes for females. 


2.1 Evidence on the Feasibility of the Identifying Assumption 


Strategy 


In the absence of a randomized treatment assignment, I need to assume that the percentage 
of women is not correlated with unobservables correlated with the outcome. As it involves 
unobservables, there is no definitive way to test this assumption. Nevertheless, I argue 
that the treatment is independent of potential outcomes and unobservables. I do this by 
analyzing how observables relate to treatment and extend this reasoning to the unobservables. 
Additionally, I present an autoregressive model to check for within major-by-branch gender 
peer composition trends. 

I first check how the treatment is correlated to some of the observable covariates in Table 
2 by regressing the covariate on the treatment using major-by-branch fixed effects. They are 
either not significant or practically zero, which is indicative that, at least in observables, the 
exogeneity assumption holds. Although certainly not conclusive, this is an indication that 


unobservables (conditional on major-by-branch fixed effects) might be perpendicular to the 


Lt 


treatment as well. 


Table 2: Balance tests 


Dependent variable: 


% of Female Peers (Treatment) 


STEM Non-STEM 
Age (X = 21.4) -0.0001** —0.00003 
(0.00004) (0.0001) 
Diagnostic score (X = 46.5) 0.00002** —0.00000 
(0.00001) (0.00002) 
Mothers’ education (X = 0.64) —0.001*** —0.001 
(0.0004) (0.001) 
Works (X = 0.59) —0.010*** —0.002 
(0.0004) (0.001) 
Has financial aid (X = 0.67) —0.001* —0.004*** 
(0.0004) (0.001) 
Night Shift (X = 0.3) —0.002*** —0.002 
(0.0004) (0.002) 
Observations 44,234 59,912 
Major-by-branch & Year Fixed Effects Yes Yes 
Note: *p<0.1; **p<0.05; ***p<0.01 


Errors clustered to the Year and Major level 


As shown in Figure 2, the percentages of female peers are similar year to year. Never- 
theless, the specification used is conditional on the major-by-branch, meaning that although 
the average level of percentage female peer might be predictable, the idiosyncratic change 
can be thought of as random. Therefore, as it is hard that students manipulate the per- 
centage of female peers they are exposed to by unilaterally switching to other majors (both 
because they do not know the exact number until they start classes when the enrollment 


process is finished and because most majors do not have a substitute that is accessible to the 
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students), assuming that the small changes within major-by-branch year-to-year are random 
is plausible. If this were not the case, we would see that the percentage of female peers in 
one year is a good predictor of the percentage of female peers in the next. To check if the id- 
iosyncratic variation in gender peer composition within programs has any trends, I estimate 


an autoregressive model of gender peer composition with major-by-branch fixed effects: 


% femaleirx = Bo + Bi % female; y_1)q + 9x + Eitk (4) 


Table 3 shows the results from estimating equation 4. In this context, peer composition 
in the previous year is not significantly correlated to gender peer composition in the current 
year, supporting the idea that within major-by-branch, the change of percentage of female 


peers is idiosyncratic. 


Table 3: Autocorrelation model for percentage of female peers 


Percentage female (t-1) 


Percentage female t —0.153 
(0.122) 
Major-by-branch & Year Fixed Effects Yes 
N 86,704 
Note: *p<0.1; *p<0.05; ***p<0.01 


STEM students sample. Errors clustered to the Year and major-by-branch level 


2.2 Relevant Variation in the Treatment for Interpretation 


To validly interpret the results of my fixed effects model, I need to find a plausible hypo- 
thetical change in the percentage of female peers supported in the data. To do so, I will use 
within-unit variation to motivate counterfactuals when discussing the substantive impact of 
the treatment (Mummolo & Peterson, 2018) 


To identify a benchmark for an increase of female percentage that is feasible and has 
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Figure 3: Within major-by-branch ranges of treatment 


support in the data, I analyze the within-unit ranges of treatment for the major-by-branch 
units. Figure 3 shows the ranges of treatment: the mean difference between the maximum 
and minimum percentage of female peers within the 106 STEM Major-by-Branch units, 
shown in blue, is 11 percentage points (p.p.), with a standard deviation of 6 p.p. To interpret 


the results, I will use 10 p.p. as a benchmark for a feasible change in treatment. 


3 Results 


Table 8 reports the estimation of the fixed effects model described in equation 3, for two out- 
comes: first-year dropout and standardized GPA. The sample includes 104,284 observations 
for dropout and 102,571 for standardized GPA, divided into STEM Major and Non-STEM 
Major students. Both regressions include student-level controls (age, diagnostic test score, 


mother’s education, working status, funding status, and night shift) and year fixed effects, 
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Major-by-Branch fixed effects, and Major-by-Branch linear time trends. The errors are clus- 
tered at the Major-by-Branch and year level. The interaction term allows comparing the 


coefficient for percentage female for males and females in each sample. 


Table 4: Estimates of the Effect of Percentage Female on Dropout 


GPA(SD) Dropout GPA(SD) Dropout 
STEM Majors Non-STEM Majors 
Perc. Female 0.534** —0.161* —0.079 0.055 
(0.250) (0.091) (0.146) (0.059) 
Male —0.075* —0.013 —0.346*** 0.048*** 
(0.041) (0.017) (0.028) (0.014) 
Perc. Female:Male —0.458** 0.115 0.066* —0.001 
(0.199) (0.078) (0.040) (0.027) 
Major-by-Branch and Year F.E. Yes Yes Yes Yes 
Major-by-Branch Time Trend Yes Yes Yes Yes 
Observations 43,356 44,234 59,085 59,912 
Note: *p<0.1; **p<0.05; ***p<0.01 


Errors clustered to the Major-by-Branch level 


Using dropout as the outcome, the coefficient of interest for students in STEM majors 
is negative and significant at the 10% level. An increase of 10 p.p. on the percentage 
of women within a STEM Major-by-Branch causes a decrease of 1.6 percentage points in 
female students’ dropout rate. This reduction represents a 9.5% decrease in dropout rates 
for females in STEM majors. 

In contrast to STEM major students, for students in Non-STEM majors, the coefficient of 
interest on dropout is small and non-significant. This suggests that gender peer composition 
does not affect dropout or GPA for Non-STEM major students. 

In the case of standardized GPA as the outcome, the estimations showed in Table ?? are 
consistent with what is observed for dropout. For STEM majors, the coefficient is positive 


and significant to the 5% level. Here, a 10 p.p. increase of the percentage of females within 
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Table 5: Estimates of the Effect of Percentage Female on Dropout and GPA (SD) 


GPA(SD) Dropout GPA(SD) Dropout 
STEM Majors Non-STEM Majors 
Perc. Female 1.293** —0.514** 0.274** —0.074** 
(0.227) (0.070) (0.097) (0.028) 
Male —0.089** —0.001 —0.301** 0.048** 
(0.030) (0.011) (0.029) (0.011) 
Perc. Female:Male —0.313* 0.112* 0.100* —0.001 
(0.146) (0.051) (0.048) (0.018) 
Mean of outcome for women —27.625 —2.770 1.220 —0.505 
Major-by-Branch F.E. Yes Yes Yes Yes 
Year F.E. No No No No 
Major-by-Branch Time Trend No No No No 
Observations 50756 52316 75332 77062 
Note: +<0.1; *<0.05; **<0.01 


Errors clustered to the Major-by-Branch level. 


the Major-by-Branch is related to a 0.05 standard deviation increase in GPA. 

Some could argue that the problem is not some particular gender dynamic in STEM 
majors but male-concentrated majors. If this were the case, we would observe a positive 
effect of a higher female percentage for women on male-concentrated non-STEM majors. 
This can be explored with the data, comparing estimates of the main specification for male 
concentrated STEM majors versus male concentrated non-STEM majors °. The results of 
estimating the main specification (equation 3) for STEM male-concentrated majors and non- 
STEM male concentrated majors are shown in Table 9 for standardized GPA and Table 10 
for dropout. 

In the estimations for STEM male-concentrated majors, we observe that the sign of the 


effect reflected in the coefficient of the percentage of female peers remains the same and is 


°Male concentrated majors are defined as majors that have less than 30% female students in their student 
body. I will run robustness checks for this benchmark in future work. 
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Table 6: Estimates of the Effect of Percentage Female on Dropout and GPA (SD) 


GPA(SD) Dropout 
STEM Majors 


GPA(SD) 


Dropout 


Non-STEM Majors 


Perc. Female 0.459* —0.177* 0.081 0.011 
(0.212) (0.070) (0.099) (0.029) 
Male —0.098** —0.012 —0.325** 0.041** 
(0.028) (0.010) (0.029) (0.011) 
Perc. Female:Male —0.358* 0.116* 0.065 0.007 
(0.139) (0.048) (0.047) (0.017) 
Mean of outcome for women —9.806 —0.954 0.361 0.075 
Major-by-Branch F.E. Yes Yes Yes Yes 
Year F.E. Yes Yes Yes Yes 
Major-by-Branch Time Trend No No No No 
Observations 50756 52316 75332 77062 
Note: +<0.1; *<0.05; **<0.01 


Errors clustered to the Major-by-Branch level. 


still significant for standardized GPA. On the other hand, for Non-STEM male-concentrated 


majors, the sign flips for both coefficients, and they are statistically insignificant. Therefore, 


having a higher percentage of female peers is associated with higher GPA, and lower dropout 


is not an artifact of majors being mostly male-concentrated but of other mechanisms that 


seem to be unique to STEM majors. 


Le 


Table 7: Estimates of the Effect of Percentage Female on Dropout and GPA (SD) for STEM 
Majors 


Dropout GPA (SD) 
(1) (2) (3) (4) (5) (6) 
Perc. Female -0.514** -0.177* -0.158* 1.293**  0.459* 0.566** 
(0.070) (0.070) (0.069) (0.227) (0.212) (0.192) 
Male -0.001 -0.012 -0.015 -0.089** -0.098** -0.089** 
(0.011) (0.010) (0.010) (0.030) (0.028) (0.027) 
Perc. Female:Male 0.112* 0.116* 0.127** -0.313* -0.358* -0.412 ** 
(0.051) (0.048) (0.048) (0.146) (0.139) (0.127) 
Mean outcome for women 0.186 0.186 0.186 -0.047 -0.047 -0.047 
Major-by-Branch F.E. Yes Yes Yes Yes Yes Yes 
Controls No Yes Yes No Yes Yes 
Major-by-Branch time trends No No Yes No No Yes 
Observations 52316 52316 52316 50756 50756 50756 
Note +<0.1; *<0.05; **<0.01. Errors clustered to the 


Major-by-Branch-by- Year level. Controls: age, diagnostic scores, mother’s education, 
working status, financial aid status, night education, and year fixed effects. 
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Table 8: Estimates of the Effect of Percentage Female on Dropout and GPA (SD) for non- 


STEM Majors 


Dropout GPA (SD) 
(1) (2) (3) (4) (5) (6) 


Perc. Female 


Male 


Perc. Female:Male 


-0.074* 0.011 0.019 0.274** 0.081 0.046 
(0.028) (0.029) (0.028) (0.097) (0.099) (0.096) 


0.048** 0.041 0.037 -0.301** -0.325** -0.311** 
(0.011) (0.011) (0.010) (0.029) (0.029) (0.028) 


-0.001 0.007 0.015 0.100* 0.065 0.040 
(0.018) (0.017) (0.017) (0.048) (0.047) (0.046) 


Mean outcome for women 
Major-by-Branch F.E. 
Controls 

Major-by-Branch time trends 
Observations 


0.146 0.146 0.146 0.225 0.225 0.225 


Yes Yes Yes Yes Yes Yes 
No Yes Yes No Yes Yes 
No No Yes No No Yes 


77062 = 77062) 77062 = 75332 75332 79332 


Note 


*+<0.1; *<0.05; **<0.01. Errors clustered to the 


Major-by-Branch-by- Year level. Controls: age, diagnostic scores, mother’s education, 


working status, financial aid status, night education, and year fixed effects. 
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Table 9: Estimates of the Effect of Percentage Female on Standardized GPA on Male- 


concentrated Majors 


Standardized GPA 
STEM Male concentrated 


Non-STEM Male concentrated 


% Female 0.514** —0.033 
(0.259) (1.267) 
Male —0.101** —0.187 
(0.043) (0.150) 
% Female:Male —0.256 —0.273 
(0.222) (0.582) 
Major-by-Branch and Year F.E. Yes Yes 
Major-by-Branch Time Trends Yes Yes 
Observations 40,372 4,309 


Note: 


SHO) 0.05)" p00) 


Errors clustered to the Year-Major-by-Branch level 


Table 10: Estimates of the Effect of Percentage Female on Dropout on Male-Concentrated 


Majors 
Dropout 
STEM Male Concentrated | Non-STEM Male Concentrated 

% Female —0.115 0.542 

(0.147) (0.459) 
Male —0.0001 0.065 

(0.019) (0.069) 
% Female:Male 0.018 —0.075 

(0.124) (0.296) 
Major-by-Branch and Year F.E. Yes Yes 
Major-by-Branch Time Trends Yes Yes 
Observations 41,182 4,386 


Note: 


*p<0.1; **p<0.05; ***p<0.01 


Errors clustered to the Major-by-Branch level 
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4 Mechanisms 


There are several potential explanations of why women improve their academic outcomes 
when surrounded by more women. One channel that has been studied extensively is the role 
model effect (Porter & Serra, 2020). If we think of teachers as role models, and if students 
identify themselves more with same-sex role models, performance may be enhanced when 
students are assigned to a same gender teacher (Dee, 2007). In this section, I explore the 
role-model channel, through which having a higher percentage of female instructors may 
interact with gender peer effects. If that was the case, then the gender peer effect I observe 
in this context is different depending on the percentage of female instructors of the student. 
If having more female instructors is a substitute for having a higher percentage of female 
peers, then the coefficient on the interaction term and the share of female coefficient would 
have opposite signs. If, on the other hand, these are complements, then both coefficients 
would have the same sign as having a higher proportion of female instructors would magnify 
the effect of having a higher share of female peers. To shed light on this idea, I estimate 
the original model, adding an interaction term between the percentage of female peers and 


percentage of female instructors: 


Yitk = Bo + Bi % female, + BoMale iz, + B3% female, * Male. + Ba% femaleinstructors jth+ 
B5% female, * % femaleinstructors tk + BoXitk + Vr + 5% + Wevearst + Eitk 


(5) 


Here, the coefficient of interest is the marginal effect of the percentage of female peers 


for women, which will include a component of the percentage of female instructors: 


OV ix 


————— = 8, + 83 *0+ 85% femaleinstructors jz, (6) 
0% female, 


And its standard errors are composed by the variance of percentage of female peers 
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(% female,,), variance of female instructors (%femaleinstructors,;), the covariance between 


both variables, and the value of the share of female instructors % femaleinstructors jx: 


6 a, = \/var(1) + % femaleinstructors?, * var(Bs) +2 * %femaleinstructorsigy * cov(B1, 63) 
Th female, 
(7) 


Given the standard error of the marginal effect shown in equation 7, if the covariance 
between two variables is negative, then it is entirely positive that the linear combination 
is significant for substantively relevant values of % femaleinstructors;,, even if the model 
parameters are insignificant ©. 

The estimates of model 5 can be found in Table 11. For easier interpretation, I only 
show the coefficients relevant to female students. When the %Femalelnstructors takes the 
value of the mean percentage of female instructors that students in STEM have (30%), 
the marginal effect of %Female in dropout is -16.8 percentage points (—0.26 + 0.3 « 0.312), 
and it is significant at the 10% level. The opposite signs suggest that as students have 
a higher percentage of female instructors, the gender peer effect caused by having more 
female peers decreases: if instead of having 30% female instructors, the students had 35% 
of female instructors, the marginal effect of %Female in dropout would drop from -16.8 to 
-15.2 percentage points (—0.26 + 0.35 *0.312). 

In the case of GPA, we also observe that the percentage of female instructors seems to be a 
substitute for having a higher share of female peers. Here, when the %FemalelInstructors takes 
the value of the mean percentage of female instructors that students in STEM have (30%), 
the marginal effects of %Female in the GPA is 0.54 standard deviations (0.645 — 0.3 « 0.344), 
and it is significant at the 5% level. The opposite signs suggest that as students have 
a higher percentage of female instructors, the gender peer effect caused by having more 


female peers decreases: if instead of having 30% female instructors, the students had 35% of 


°For more information and practical examples of this, see (Brambor, Clark, & Golder, 2006) 
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female instructors, the marginal effect of %Female in GPA would decrease from 0.54 to 0.52 


percentage points (0.645 — 0.35 « 0.344). 


Table 11: Percentage of Female Instructor Interaction 


Dropout GPA (Standardized) 
New Model Original Model New Model Original Model 
% Female —0.261** —0.161* 0.645* 0.534** 
(0.119) (0.091) (0.336) (0.250) 
% Female Instructors —0.058 0.068 
(0.045) (0.121) 
% Female:% Female Instructors 0.312* —0.344 
(0.175) (0.499) 
Major-by-Branch & Year F.E. Yes Yes Yes Yes 
Major-by-Branch Time Trends Yes Yes Yes Yes 
N 43,645 44,234 42,782 43,356 


5 Discussion 


Note: *p<0.1; *p<0.05; ***p<0.01 
Errors clustered at the Major-by-Branch and Year level 


The results presented in this paper support the hypothesis that there are gender peer effects 


in the vocational education system. In particular, having a higher percentage of female peers 


positively affects students in STEM majors, decreasing women’s dropout rates and GPA in 


vocational post-secondary education. Furthermore, a higher percentage of female peers on 


men decreases dropout, but on a smaller scale. The evidence presented in this paper suggests 


that both women and men benefit from having a higher percentage of female peers. This has 


a two-fold implication: first, it supports the scholarship that has established that a higher 


share of female peers is associated with better outcomes (Lavy & Schlosser, 2011). Second, 


it confirms Hoxby (2001)’s finding that both males and females benefit from having a higher 
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percentage of female peers, challenging recent evidence that students’ outcomes are hindered 
by having a higher share of opposite gender schoolmates (Hill, 2017). 

In terms of policymaking, this paper presents evidence that gender peer effects exist 
in vocational education and that they might point to an intervention path to “stop the 
leaking” in the STEM sector. The mechanism analysis also suggests that increasing the 
percentage of female instructors could make female students’ outcomes better by substituting 
the percentage of female peers. Although post-secondary vocational institutions cannot 
directly control the gender composition of cohorts in this context, they do have discretion in 
instructor selection. Therefore, these results suggest that increasing the percentage of female 
instructors is an avenue to improve female students’ outcomes when increasing the percentage 
of female peers is not feasible. Nevertheless, more causal evidence on this mechanism is 
necessary to affirm that role models will benefit women in this context. As scholars and 
institutions develop analyses of gender peer composition, collecting data and using causal 
inference methodologies will allow identifying levers to avoid gender polarization in STEM. 
From a theory perspective, this paper proposes that in the case of vocational education, 
having a higher percentage of female students represents a Pareto improvement: both women 
and men benefit from it, or at the very least are not harmed by it. Although the literature 
has studied the effect on women extensively, there is little evidence on the effect that gender 
composition has on men, an important point when thinking about the general welfare of 
students. 

A big question when implementing policies that improve women’s outcomes in education 
is if there is a trade-off between improvements for women and men. The estimations presented 
in this paper provide evidence that for STEM majors, an increase of females within major- 
by-branches would benefit women and would not harm men. This is an essential point for 
policy-makers and institutions that want to push for policies and strategies to increase female 
participation in STEM programs. 


This paper builds on the gender peer effect literature, providing evidence for a novel 
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context (a middle-income country in Latin America) for an education sector that has not been 
thoroughly studied: post-secondary vocational education. It provides strong evidence that 
the international trends of female student achievement and peer effects hold in this context 
and that actions geared toward improving gender balance within majors can positively affect 
students’ outcomes. To understand better the actions that could improve gender balance in 
this context, the next step should be to explore the mechanisms at play that create these 


dynamics experimentally. 
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Appendices 


A 


Female Percentage 


Identifying variation plots 


Figure 4: Within major-by-branch Yearly percentage female variation 
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B Mechanism Full Estimation 


Table 12: Percentage of Female Instructor Interaction 


Table 12: Percentage of Female Instructor Interaction 


Dropout GPA (Standardized) 
New Model Original Model New Model — Original Model 

% Female —0.261** —0.161* 0.645* 0.534** 

(0.119) (0.091) (0.336) (0.250) 
% Female Instructors —0.058 0.068 

(0.045) (0.121) 
Male —0.013 —0.013 —0.077 —0.075* 

(0.018) (0.017) (0.047) (0.041) 
% Female:% Female Instructors 0.312* —0.344 

(0.175) (0.499) 
% Female:Male 0.116 0.115 —0.450** —0.458** 

(0.081) (0.078) (0.209) (0.199) 
Major-by-Branch & Year F.E. Yes Yes Yes Yes 
Major-by-Branch Time Trends Yes Yes Yes Yes 
N 43,645 44,234 42,782 43,356 


Note: *p<0.1; “p<0.05; ***p<0.01 
Errors clustered at the Major-by-Branch and Year level 
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Table 6: Estimates of the Effect of Percentage Female on Dropout and GPA (SD) 


GPA(SD) Dropout GPA(SD) Dropout 


STEM Majors Non-STEM Majors 

Perc. Female 0.459* —0.177* 0.081 0.011 

(0.212) (0.070) (0.099) (0.029) 
Male —0.098** —0.012 —0.325** 0.041** 

(0.028) (0.010) (0.029) (0.011) 
Perc. Female:Male —0.358* 0.116* 0.065 0.007 

(0.139) (0.048) (0.047) (0.017) 
Mean of outcome for women —9.806 —0.954 0.361 0.075 
Major-by-Branch F.E. Yes Yes Yes Yes 
Year F.E. Yes Yes Yes Yes 
Major-by-Branch Time Trend No No No No 
Observations 50756 52316 75332 77062 
Note: +<0.1; *<0.05; **<0.01 


Errors clustered to the Major-by-Branch level. 


dl 


Table 7: Estimates of the Effect of Percentage Female on Dropout and GPA (SD) for STEM 
Majors 


Dropout GPA (SD) 
(1) (2) (3) (4) (5) (6) 
Perc. Female -0.514** -0.177* -0.158* 1.293**  0.459* 0.566** 
(0.070) (0.070) (0.069) (0.227) (0.212) (0.192) 
Male -0.001 -0.012 -0.015 -0.089** -0.098** -0.089** 
(0.011) (0.010) (0.010) (0.030) (0.028) (0.027) 
Perc. Female:Male 0.112* 0.116* 0.127** -0.313* -0.358* -0.412 ** 
(0.051) (0.048) (0.048) (0.146) (0.139) (0.127) 
Mean outcome for women 0.186 0.186 0.186 -0.047 -0.047 -0.047 
Major-by-Branch F.E. Yes Yes Yes Yes Yes Yes 
Controls No Yes Yes No Yes Yes 
Major-by-Branch time trends No No Yes No No Yes 
Observations 52316 52316 52316 50756 50756 50756 
Note *+<0.1; *<0.05; **<0.01. Errors clustered to the 


Major-by-Branch-by- Year level. Controls: age, diagnostic scores, mother’s education, 
working status, financial aid status, night education, and year fixed effects. 
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Table 8: Estimates of the Effect of Percentage Female on Dropout and GPA (SD) for non- 
STEM Majors 


Dropout GPA (SD) 
(1) (2) (3) (4) (5) (6) 
Perc. Female -0.074** 0.011 0.019 0.274** 0.081 0.046 
(0.028) (0.029) (0.028) (0.097) (0.099) (0.096) 


Male 0.048" 0.041** 0.037 -0.301** -0.325**  -0.311** 
(0.011) (0.011) (0.010) (0.029) (0.029) (0.028) 


Perc. Female:Male -0.001 0.007 0.015 0.100* 0.065 0.040 
(0.018) (0.017) (0.017) (0.048) (0.047) (0.046) 
Mean outcome for women 0.146 0.146 0.146 0.225 0.225 0.225 
Major-by-Branch F.E. Yes Yes Yes Yes Yes Yes 
Controls No Yes Yes No Yes Yes 
Major-by-Branch time trends No No Yes No No Yes 
Observations 77062. =©77062) =677062 =——75332 (0832 75332 
Note +<0.1; *<0.05; **<0.01. Errors clustered to the 
Major-by-Branch-by- Year level. Controls: age, diagnostic scores, mother’s education, 


working status, financial aid status, night education, and year fixed effects. 
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Table 9: Estimates of the Effect of Percentage Female on Standardized GPA on Male- 
concentrated Majors 


Standardized GPA 
STEM Male concentrated Non-STEM Male concentrated 


% Female 0.514** —0.033 

(0.259) (1.267) 
Male —0.101** —0.187 

(0.043) (0.150) 
% Female:Male —0.256 —0.273 

(0.222) (0.582) 
Major-by-Branch and Year F.E. Yes Yes 
Major-by-Branch Time Trends Yes Yes 
Observations 40,372 4,309 
Note: *p<0.1; *p<0.05; ***p<0.01 


Errors clustered to the Year-Major-by-Branch level 
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Table 10: Estimates of the Effect of Percentage Female on Dropout on Male-Concentrated 
Majors 


Dropout 
STEM Male Concentrated | Non-STEM Male Concentrated 
% Female —0.115 0.542 
(0.147) (0.459) 
Male —0.0001 0.065 
(0.019) (0.069) 
% Female:Male 0.018 —0.075 
(0.124) (0.296) 
Major-by-Branch and Year F.E. Yes Yes 
Major-by-Branch Time Trends Yes Yes 
Observations 41,182 4,386 
Note: *p<0.1; *p<0.05; ***p<0.01 


Errors clustered to the Major-by-Branch level 
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Table 11: Percentage of Female Instructor Interaction 


Dropout GPA (Standardized) 
New Model Original Model New Model Original Model 
% Female —0.261** —0.161* 0.645* 0.534** 
(0.119) (0.091) (0.336) (0.250) 
% Female Instructors —0.058 0.068 
(0.045) (0.121) 
% Female:% Female Instructors 0:312* —0.344 
(0.175) (0.499) 
Major-by-Branch & Year F.E. Yes Yes Yes Yes 
Major-by-Branch Time Trends Yes Yes Yes Yes 
N 43,645 44,234 42,782 43,356 
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Note: *p<0.1; **p<0.05; ***p<0.01 
Errors clustered at the Major-by-Branch and Year level 


Table 12: Percentage of Female Instructor Interaction 


Dropout GPA (Standardized) 
New Model Original Model New Model Original Model 

% Female —0.261** —0.161* 0.645* 0.534** 

(0.119) (0.091) (0.336) (0.250) 
% Female Instructors —0.058 0.068 

(0.045) (0.121) 
Male —0.013 —0.013 —0.077 —0.075* 

(0.018) (0.017) (0.047) (0.041) 
% Female:% Female Instructors 07312" —0.344 

(0.175) (0.499) 
% Female:Male 0.116 0.115 —0.450** —0.458** 

(0.081) (0.078) (0.209) (0.199) 
Major-by-Branch & Year F.E. Yes Yes Yes Yes 
Major-by-Branch Time Trends Yes Yes Yes Yes 
N 43,645 44,234 42,782 43,356 


Note: *p<0.1; **p<0.05; ***p<0.01 
Errors clustered at the Major-by-Branch and Year level 
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Table 12: Percentage of Female Instructor Interaction 


Dropout GPA (Standardized) 
New Model Original Model New Model Original Model 

% Female —0.261** —0.161* 0.645* 0.534** 

(0.119) (0.091) (0.336) (0.250) 
% Female Instructors —0.058 0.068 

(0.045) (0.121) 
Male —0.013 —0.013 —0.077 —0.075* 

(0.018) (0.017) (0.047) (0.041) 
% Female:% Female Instructors 07312" —0.344 

(0.175) (0.499) 
% Female:Male 0.116 0.115 —0.450** —0.458** 

(0.081) (0.078) (0.209) (0.199) 
Major-by-Branch & Year F.E. Yes Yes Yes Yes 
Major-by-Branch Time Trends Yes Yes Yes Yes 
N 43,645 44,234 42,782 43,356 


Note: *p<0.1; **p<0.05; ***p<0.01 
Errors clustered at the Major-by-Branch and Year level 
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